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Abstract 

In order to elucidate the emission properties of ultraluminous X-ray sources (ULXs) during their power-law (PL) 
state, we examined long-term X-ray spectral data of IC 342 X-1 during its PL state by using our own Suzaku data 
and the archival data by XMM-Newton, Chandra, and Swift observations. The PL state of this source seems to be 
classified into two sub-states in terms of the X-ray luminosities in 0.5-10 keV: the low luminosity PL state with 
(4-6)xlG'^^ erg s~^ and the high luminosity one with ( 1.1-1. 4) x 10^° erg s^^. During the Suzaku observations 
which were made in 2010 August and 201 1 March, X-1 stayed in the low luminosity PL state. The observed X-ray 
luminosity (4.9-5.6 x 10^^ erg s^^) and the spectral shape (photon index = 1.67-1.83) slightly changed between the 
two observations. Using the Suzaku PIN detector, we for the first time confirmed a PL tail extending up to at least 
20 keV with no signatures of a high-energy turnover in both of the Suzaku observations. In contrast, a turnover at 
about 6 keV was observed during the high luminosity PL state in 2004 and 2005 with XMM-Newton. Importantly, 
photon indices are similar between the two PL states and so is the Compton j/-parameters of y ^ 1, which indicates 
a similar energy balance (between the corona and the accretion disk) holding in the two PL states despite different 
electron temperatures. From spectral similarities with recent studies about other ULXs and the Galactic black hole 
binary GRS 1915+105, IC342 X-1 is also likely to be in a state with a supercritical accretion rate, although more 
sensitive higher energy observations would be necessary to conclude. 
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1. Introduction 

Ultraluminous X-ray sources (ULXs; Makishima et al. 
2000) continue to be enigmatic objects since their discovery in 
the 1980's; see Fabbiano (2004) and Feng & Soria (2011) for 
reviews. They are bright, variable, and off-nuclear point-Uke 
X-ray sources with luminosities of >10'^^~'*^ erg s~^, which 
exceed the Eddington luminosity (iEdd) of Galactic black 
holes (BHs). Despite intensive multi-wavelength studies over 
three decades, their central engine remains an open question. 
ULXs are very likely to contain a BH, but not a super-massive 
BH, since their X-ray properties are similar to those of Galactic 
BH binaries (BHBs) in many respects; e.g., they show short- 
term flux variations and have records of exhibiting at least two 
spectral states characterized by power-law (PL) shaped spectra 
and convex-shaped (thermal) spectra (Makishima 2007; Soria 
201 1). Some ULXs have even exhibited spectral transitions be- 
tween the two (La Parola et al. 2001; Kubota et al. 2001; Isobe 
et al. 2009). These similarities indicate the same mechanism 
(i.e., gas accretion onto BHs) being responsible for electromag- 
netic radiation in ULXs. 

Some authors (e.g., Kubota et al. 2002; Mizuno et al. 
2007; Makishima 2007; Isobe et al. 2009) claimed that the PL 
and the convex spectral states of ULXs may correspond to the 
very high state (VHS; Miyamoto et al. 1991) and an extremely 



high luminous disk state (Kubota & Makishima 2004; Isobe et 
al. 2012) described by a slim disk model (Abramowicz et al. 
1988) of Galactic BHBs, respectively. This interpretation is 
based on the following important observational facts: (i) Some 
observed properties of ULXs are consistent with the theoretical 
study based on the slim disk model (Watarai et al. 2000). For 
example, the derived temperature profile of the disk is flatter 
than that of the standard disk (Vierdayanti et al. 2006; Tsunoda 
et al. 2006), and the apparent innermost disk radius is inversely 
proportional to the innermost disk temperature (Mizuno et al. 
2001). (ii) The luminosity difference between the two states of 
ULXs are consistent with that between the slim disk state and 
the VHS of Galactic BHBs. Such ULXs were reported to be 
brighter in the convex state than in the PL state (La Parola et 
al. 2001; Kubota et al. 2001; Isobe et al. 2009). 

However, some concerns remain. First, the PL state of 
ULXs shows a wider range of the photon indices; for example, 
r = 1.0-1.5 (M51 source-69 by Yoshida et al. 2010), 1.6-2.0 
(NGC 1365 X-1 by Soria et al. 2009), and 2.0-2.7 (NGC 5204 
X-1 by Kajava & Poutanen 2009). Such various photon in- 
dices have not been reported in the VHS of Galactic BHBs. For 
this concern, we are unable to rule out the possibility that the 
energy range of typical X-ray instruments (e.g., 0.5-10 keV) 
corresponds to the different portion of the intrinsic X-ray spec- 
trum in various types of ULXs (and BHBs), and that it makes 
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the variety in the photon index of ULXs. The second con- 
cern is that a spectral turnover at 3.5-7 keV (e.g., Kubota et 
al. 2002; Mizuno et al. 2007; Gladstone et al. 2009), which is 
often accompanied with the PL spectra of ULXs, is found at 
much lower energy than that of Galactic BHBs in the VHS (at 
about several tens keV; Kubota & Done 2004). These differ- 
ences seem to imply that the PL state of ULXs and the VHS of 
Galactic BHBs could be distinct in nature. 

Obviously, investigating the characteristics of ULXs in vari- 
ous PL states is a key to unveiling the origin of ULXs, though 
they have not been so extensively discussed so far, compared 
with those in the convex state. The suitable target to achieving 
our purpose is the unique ULX X-1 in the nearby spiral galaxy 
IC 342 (figure la) at 3.3 Mpc (Saha et al. 2002). While this is 
the first ULX exhibiting the transition from the PL state to the 
convex one (Kubota et al. 2001; Kubota et al. 2002), it stayed 
in the PL state in most of the previous observations (Roberts et 
al. 2004; Feng & Kaaret 2009; Mak et al. 2011). Further two 
seemingly distinct sub-states have been reported; the high lu- 
minosity PL state and the low luminosity PL state, where the 
former is brighter than the latter by a factor of 2 to 3 (Feng & 
Kaaret 2009; Mak et al. 2011). 

In order to elucidate the nature of the PL state of the source, 
we discuss the PL spectral states and the long-term spectral 
variabilities of IC342 X-1 assembling all the recent available 
data sets, which include our new data taken with Suzaku and 
unpublished data taken with Swift. A particular attention is 
paid to the Suzaku data, since it has a detector sensitive above 
10 keV (see § 2.1 for details), thereby being able to determine 
the existence of a turnover at a hard X-ray band in the PL state. 

This paper is composed of the following sections: In § 2, 
we summarize the observations and data reduction. In § 3 and 
§ 4, we describe the results of timing and spectral analysis, 
respectively. We develop discussion for the PL states of IC 342 
X-1 in § 5, and conclude this paper in § 6. 

2. Observations and Data Reduction 

Table 1 shows the assembled data sets. We labeled the ob- 
servations as Sul-2 for Suzaku, XMl^ for XMM-Newton, 
Chl-2 for Chandra, and Swl-3 for Swift. The list includes 
all the available data sets except for a Chandra observation 
(ObsID = 7069), because the target source suffers a severe pile- 
up with a pile-up fraction of more than 20%. Throughout this 
paper, we use events in the 14-20 keV band for Suzaku PIN, 
and in 0.5-10 keV for the other instruments. 

2.1. Suzaku 

The Suzaku observatory (Mitsuda et al. 2007) has two 
operating science instruments. One is the X-ray Imaging 
Spectrometer (XIS; Koyama et al. 2007) covering an energy 
range of 0.2-12 keV, while the other is the Hard X-ray Detector 
(HXD; Takahashi et al. 2007; Kokubun et al. 2007) covering 
10-600 keV. 

We observed IC342 with Suzaku twice on 2010 August 
(Sul) and 2011 March (Su2) for about 70 ks each (table 1). 
We aimed at the middle point of X-1 and the other ULX (X-2) 
in order to observe the two ULXs within the XIS field of view 



Table 1. Observation log. 



Data 


Observatory 


ObsID 


Date 


^cxp 


label 








(ks) 


Sul 


Suzaku 


705009010 


2010-08-07 


74.4/64.1 


Su2 




705009020 


2011-03-20 


68.5/74.3 


XMl 


XMM-Newton 


0093640901 


2001-02-11 


9.7/ 5.1 


XM2 




0206890101 


2004-02-20 


23.4/17.8 


XM3 




0206890201 


2004-08-17 


23.4/17.1 


XM4 




0206890401 


2005-02-10 


18.6/15.1 


Chi 


Chandra 


2916 


2002-04-29 


9.3 


Ch2 




2917 


2002-08-26 


9.9 


Swl 


Swift 


00035474001 


2006-01-07 


7.2 


Sw2 




00031987001 


2011-06-14 


3.0 


Sw3 




00031987002 


2011-06-15 


9.4 



* XIS (left) and HXD PIN (right) exposures for Suzaku, KPIC MOS 
(left) and EPIC pn (right) exposures for XMM-Newton, ACIS exposure for 
Chandra, and XRT exposure for Swift. 



(FoV). The data in these observations were taken with the nor- 
mal clocking mode with a frame time of 8 s for XIS, and with 
the normal mode with the time resolution of 61 /is for HXD. 

In order to create new cleaned event files, we reprocessed 
the data using the software package HEASoft' version 6.10 
with the X-ray telescope calibration database (CALDB) ver- 
sion 20100730, the XIS CALDB version 20110210, and the 
HXD CALDB version 20101202. The data were screened out 
during the South Atlantic Anomaly passages (SAA), within 
436 s after exiting from the SAA, with the Earth elevation angle 
of below 5°, with the day-Earth elevation angle of below 20° 
(for XIS), or at low cutoff rigidity regions (<6 GV for HXD). 

2.1.1. XIS 

XIS is composed of three front-illuminated (FI) and one 
back-illuminated (BI) CCD devices with a FoV of 17f8x 17^8. 
Since one of the three FI devices became dysfunctional due 
to a putative micrometeorite hit in 2006 November, we used 
the remaining two FI and one BI devices for our analysis. We 
merged the two FI data because the devices have nearly identi- 
cal responses. 

Figures 1 (b) and (c) show Suzaku XIS three-color images of 
IC 342. We can see the three bright sources (X-1, X-2, and X-3 
named by Okada et al. 1998) and two relatively faint sources 
(sources-6 and 11; Bauer et al. 2003). X-2 and X-3 are re- 
spectively harder and softer than X-1 in the XIS range. Since 
hard sources may contribute to the PIN spectrum (see § 2.1.2), 
we also investigated the properties of X-2 in the present pa- 
per In the equinox J2000.0, the peak positions of X-1 and X-2 
were respectively determined at (R.A., Dec.) = (03'^45™55*', 
+68°04'59") and (03''46"'15", +68°11'17"), which are con- 
sistent with previous studies (Kubota et al. 2001; Kong 2003; 
Bauer et al. 2003; Mak et al. 201 1). 

Source events were extracted from a circle of an adaptively 
chosen radius of 130", 130", 140", and 120", for X-1 in Sul, 
X-1 in Su2, X-2 in Sul, and X-2 in Su2, respectively (fig- 
ures lb and c). Whereas, background events were extracted 
from a common region. Here, we masked circles with a radius 
of 1' around the two relatively faint sources-6 and 11. Other 
two relatively faint sources near each ULX (source-9 for X- 
1 and source-12 for X-2; Kong 2003; Bauer et al. 2003) are 

' See http://heasarc.gsfc.nasa.gov/docs/software/lheasoft/ for details. 
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Fig. 1. Optical and X-ray images of IC 342. (a) Second Digitized Sky Survey three-color image. Red, green, and blue represent the 850, 650, and 480 nm 
band intensities, respectively. Pluses show the positions of the two ULXs. (b)-(c) Smoothed Suzaku XIS thi'ee-color images obtained in (b) 2010 August 
and (c) 201 1 March. Red, green, and blue represent the 0.5-2, 2—4, and 4-10 keV bands, respectively. Solid and dashed circles indicate the source and 
background extraction regions, respectively. The other numbers are each source name (X-3, sources-6, and 11). Solid squares show the FoV of XIS. 



not resolved in the XIS images, but make no signification con- 
tamination. Source-9 has 10% of the X-ray flux of X-1, while 
source-12has 3% of X-2. For the spectral analysis, we selected 
only the events with the standard ASCA grade set of 0, 2, 3, 4, 
and 6. 

We calculated their redistribution matrix functions (RMFs) 
using the xisrmf gen tool. Whereas, the ancillary response 
files (ARFs) were simulated by the xissimarfgen tool 
(Ishisaki et al. 2007). The influence by relatively faint sources 
was removed. 

2.1.2. HXDPIN 

HXD consists of two components: PIN and GSO respec- 
tively sensitive to 10-70 keV and 40-600 keV X-rays. In this 
paper, we focus on the PIN detector only, because the object is 
too faint for a GSO detection. PIN is a non-imaging detector, 
which has a full width at zero intensity (FWZI) view of ^70' 
square. The effective area monotonically decreases as the dis- 
tance increases from the field center. 

If PIN has a high temperature, the data are affected by in- 
creased noise. In order to deal with this situation, we further 
screened PIN events by applying the temperature threshold for 
observations in the epoch 10 (2010 December 1 to 2011 May 
24)-. The threshold excludes a noisy energy band of each PIN 
unit. We used the detector response files distributed by the in- 
strument team. 

2.2. XMM-Newton 

The details of the instruments and the observations can be 
found in Feng & Kaaret (2009), Kajava & Poutanen (2009), 
and Mak et al. (2011). We used the Science Analysis System 
(SAS) version 11.0.0 (Gabriel et al. 2004) for reprocessing the 
data, extracting events, and generating response files. We used 
the calibration files latest as of 2011 April. We excluded the 
intervals with a high background rate, which is defined as the 
10-15 keV count rate of the entire array larger than the average 

^ See http://heasarc.nasa.gov/docs/suzaku/analysis/pinepochs.html for de- 
tails. 



by more than 3a. The events in the energy range are dominated 
by the background (Read & Ponman 2003). 

A common background region devoid of bright X-ray 
sources was adopted. We extracted source events from a circle 
of adaptively chosenradius of 45" forXMl, 55" forXM2, 55" 
for XM3, and 65" for XM4. For the spectral analysis, we se- 
lected the MOS events with PATTERN < 12, #X]V1MEA_EM, 
and FLAG = 0, and the pn events with PATTERN < 4, 
#XMMEAXP, and FLAG = 0. The RMF and ARE were gen- 
erated with the SAS tools rmf gen and arf gen, respectively. 

2.3. Chandra 

The details of the instrument and the observations can be 
found in Roberts et al. (2004). We reduced the data using 
the Chandra Interactive Analysis of Observations (CIAO) ver- 
sion 4.3 with the CALDB version 4.4.2. We then extracted the 
events (with grade 0, 2, 3, 4, and 6), and constructed the energy 
spectra using the ACIS Extract package (Broos et al. 2002) ver- 
sion 2009-12-01. The source events were accumulated from 
a region around each source encircling 90% of photons of a 
point-like source. The background events were from an annu- 
lus around each source. The RMF and ARF were generated 
with the CIAO tools mkrmf and mkarf , respectively. 

2.4. Swift 

The Swift observatory (Gehrels et al. 2004) is equipped with 
some instruments. One of them is the X-ray telescope (XRT; 
Burrows et al. 2005) with an X-ray CCD device sensitive at 
0.2-10 keV. The FoV (23f6x23.'6) covers the whole of IC 342. 
In the observations, the photon counting mode with a time res- 
olution of 2.5 s was used. 

We processed the data sets based on the CALDB ver- 
sion 20110513 using HEASoft version 6.10. In order to in- 
crease the photon statistics, two data sets (Sw2 and Sw3) with 
very close observation dates were merged. Hereafter we call 
the merged observation Sw2+3. Source and background re- 
gions were determined by the method same as for XMM- 
Newton, and we then extracted the spectra from grade 0-12. 
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We used the RMF file in the CALDB, and we created the ARF 
using the xrtmkarf tool. 

3. Timing Analysis and Results 

Figure 2 shows the background-subtracted Suzaku XIS (FI) 
light curves and the hardness ratio (the 3.0-10.0 keV band 
count rate to the 0.5-3.0 keV band one) variations for each 
ULX. No significant short-term variation was found in the light 
curves nor the hardness ratio variations at 99.9% significance, 
indicating that the spectra did not change within the individual 
exposure. We thus stacked all the events in each observation to 
construct spectra. 

For X-1, the time-averaged count rate and hardness ratio 
changed at most by several percent between the two Suzaku 
observations. The time-averaged count rate was measured to 
be 0.083±0.001 counts s^^ (Sul) and0.086±0.001 counts s^^ 
(Su2), while the time-averaged harness ratio was 0.66±0.02 
(Sul) and 0.71 ±0.02 (Su2). For X-2, the count rate 
became to half, from 0.203±0.002 counts s^^ (Sul) to 
0.093±0.001 counts s^^ (Su2). Its hardness ratio (1.63±0.03 
for Sul and 1.74±0.05 for Su2) also changed only by several 
percent like X- 1 , but was clearly higher than X- 1 as mentioned 
from the three-color images (figure 1). Here, the errors are the 
la values. 

The other data sets with XMM-Newton, Chandra, and Swift 
exhibit no short-term variation in each exposure. We thus 
stacked all the events in each observation to construct spectra. 

4. Spectral Analysis and Results 

First, we explain spectral fitting models in § 4. 1 . Next, we 
report our Suzaku results based on the following three steps: 
In § 4.2.1, we fit the spectral models to the XIS spectra of X- 
1 and X-2, and describe the results. Note that the spectra of 
X-2 are analyzed in order to examine the contribution to the 
PIN spectrum. In § 4.2.2, we estimate the source signals above 
10 keV with the PIN detector by inspecting the NXB accu- 
racy. In § 4.2.3, the XIS and PIN spectra are jointly analyzed. 
Finally, we present the fitting results of the other satellites in 
§4.3. 

4.1. Fitting Models 

We used the X-ray spectral fitting package XSPEC ver- 
sion 12.6.0 for the spectral analysis. We applied two interstellar 
extinctions: One is the Galactic extinction represented by the 
wabs model (Morrison & McCammon 1983) with a hydro- 
gen column density fixed at the Galactic value toward IC 342; 
Nu = 3xl02i cm-2 (Dickey & Lockman 1990; Stark et al. 
1992). The other is a thawed additional extinction for each 
ULX (as another wabs). 

All the spectra are featureless, so we first applied two sim- 
ple models: the PL model and the multi-color disk black body 
model (MCD; diskbb in XSPEC; Mitsudaet al. 1984), which 
represent PL-shaped spectra and convex-shaped spectra, re- 
spectively. These models are often employed for Galactic 
BHBs and ULXs, including the IC 342 ULXs (e.g., Kubota et 
al. 2001; Roberts et al. 2004; Mak et al. 201 1). 

If a spectrum is represented by the PL model, we further 
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Fig. 2. Background-subtracted XIS FI light curves (black) and hard- 
ness ratio variations (gray) of the two ULXs in 0.5-10 keV. The curves 
are binned with 2000 s bin~^. The hardness is calculated as the count 
rate ratio of 3-10 keV against 0.5-3 keV. 

applied two other models in order to understand the physical 
condition of the PL state. One is the cutoff-PL model, which 
approximates the Compton emission by constraining two phe- 
nomenological parameters: the photon index and the cutoff 
(turnover) energy. The improvement from the PL model is 
checked by the i^-test. We consider that the improvement 
with the >99.9% confidence level is significant. The other is a 
Comptonization model (compTT in XSPEC; Titarchuk 1994), 
which can constrain the electron temperature (kTc) and the op- 
tical depth for electrons (r) as model parameters. If the spec- 
trum has no turnover, both parameters are strongly coupled and 
cannot be constrained. 

4.2. Suzaku Results 

4.2.1. XIS spectra 

We fitted the FI and BI spectra simultaneously. The energy 
range of 1.8-2.0 keV is excluded due to calibration uncertain- 
ties. In addition, for X-2, we did not use data below 1.5 keV 
to characterize the hard X-ray spectrum, since a soft excess re- 
ported by Feng & Kaaret (2009) and Mak et al. (2011) was 
recognized in the Suzaku data. 

The Suzaku XIS spectra and the best-fit models are shown in 
figure 3, and all the best-fit parameters (90% confidence level) 
are summarized in table 2. The parameters iV^, F, kT^, and 
i^cut indicate the absorption column density additional to that 
in our Galaxy, the photon index, the innermost disk tempera- 
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Table 2. Best-fit parameters for the Suzaku XIS spectra of the two ULXs. 

Name Data Model 1^ F or t -Ecut, kTc, or fcTin Luminosity* Red-x^ 

label (10^^ cm~^) (keV) (10='^ erg s~^) (d.o.f.) 

Sul PL 0.36 ±0.04 L83±0.04 ~ 5.6±0.1 1.16(183) 

cutoff-PL 0.27 ±0.06 1.5 ±0.2 lltf S.llJ^J;^ 1.11(182) 

compTT 0.32+°°^ 7.1"^ 2.4+ S.Slo'g 1.12(181) 

MCD <0.02 ••• L85±0.05 4.50 ±0.07 1.70(183) 

4.9±0.1 1.15(175) 

>17 4.8t[]:i 1.16(174) 

2.8'f 4.8±0.1 1.17(173) 

1.94 ±0.05 4.33 ±0.07 2.25(175) 

X-2 Sul PL 3.6 ±0.1 2.49 ±0.04 ■•■ 34 ±2 1.50(396) 

MCD 1.82 ±0.08 — ^ 1.78 ±0.03 15.0 ± 2 1.05(396) 

~Su2 PL 2.2 ±0.2 1.71^0 06 ~ lO.ltot 1.19(175) 

MCD 1.1±0.1 .^-"^^Q-.^ .^-^^^^-^ 

* the 0.5-10 keV luminosity for the PL, the cutott-PL, and the Comptonization models, and the bolometric luminosity tor the MCD model, 
t Only the best-fit value is described for r and kTe, because both parameters are strongly coupled. 



Su2 PL 0.22 ±0.04 1.67 ±0.04 

cutoff-PL Q.2ltli% L63l°°« 

compTT 0.22 ± 0.02 6.9''^ 
MCD < 0.003 
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Fig. 3. Suzaku spectra (pluses) and the best-fit models (solid histogram) of (a-b) X-1 and (c-d) X-2. Bottom panels show the residuals from the best-fit 
by a PL and MCD. 
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ture, and the turnover energy, respectively. The 0.5-10 keV lu- 
minosity Lx for the PL, the cutoff-PL, and the Comptonization 
models is calculated as AnD^f^. Assuming the disk inclina- 
tion angle of 60°, the bolometric luminosity Lboi for the MCD 
model is calculated as 27rZ?^/^^j(cos 60°)~^. Here, and 
are respectively the absorption-corrected flux in the 0.5- 
10 keV and the 0.01-100 keV ranges, and D is the distance to 
IC 342 (3.3 Mpc). 

The X-1 spectra were unaccountable by the MCD model 
(reduced-x^ = 1.70 for Sul and 2.25 for Su2), whereas 
they were successfully fitted by the PL model (reduced- 
= 1.15-1.16). Based on the derived parameters (ta- 
ble 2), the spectral shape of X-1 slightly changed be- 
tween the two observations; in the Sul observation, the 
X-ray luminosity and the photon index were respectively 
Lx = (5.6±0.1)xl039 erg s'^ and r = 1.83±0.04, while 
in Su2, they were Lx = (4.9±0.1)x 10^^ erg s^^ and 
F = 1.67±0.04. These parameters are within the typical values 
(F = 1.6-2.0) measured previously from this ULX; see, e.g., 
Kubota et al. (2001) with ASCA, Feng & Kaaret (2009) with 
XMM-Newton, and Roberts et al. (2004) with Chandra. 

The cutoff-PL model provided a fit goodness similar to that 
by the PL model, but we conclude that both X-1 spectra in Sul 
and Su2 do not have a turnover at least below 10 keV, because 
the turnover energy (E'cut = H-^^ keV for Sul and >17 keV 
for Su2; table 2) is determined mostly above the XIS fitting 
range (0.5-10 keV). In the following (§ 4.2.3), we confirm this 
conclusion by using the PIN data. 

The Comptonization model (compTT) fitting is also consis- 
tent with the result by the cutoff-PL model (i.e., no turnover), 
because kTc and r were strongly coupled in both spectra, 
shown in § 5.1.2 for details. The seed photon temperature was 
fcTscod < 0.4 keV in both spectra. 

Unlike X-1, the X-2 spectrum in Sul cannot be explained 
by the PL model (reduced-x^ = 1.50). Instead, it has a con- 
vex shape represented by the MCD model with a reduced-x^ 
of 1.05 (figure 3c). On the other hand, the X-2 spectrum in 
Su2 was reproduced by both models, although the MCD model 
provided a slightly good fitness (figure 3d; reduced-x^ =1.19 
by PL and 1.08 by MCD). The source spectrum dramatically 
changed between the two observations; the luminosity halved 
from Lboi = (1.50±0.02)xl04° to (8.9±0.2)x lO^^ erg s-\ 
and its disk temperature increased (fcTin = 1.78±0.03 to 
2.7±0.1 keV). Such a spectral behavior is opposite to that of 
Galactic BHBs in the canonical high-soft or slim disk states. 
Its physical interpretation is put aside for the forth-coming pa- 
per 

4.2.2. Significance of the source signal with PIN 

For the detailed PIN spectral analysis, we need to consider 
two background sources: the non-X-ray background (NXB) 
and the cosmic X-ray background (CXB). We adopted NXB 
spectra provided by the instrument team (Fukazawa et al. 
2009). Meanwhile, for CXB, we simulated spectra based on 
the HEAO-1 model (Boldt 1987). The count rates of the ob- 
served PIN signal, the modeled NXB, the simulated CXB, and 
the source signal are summarized in table 3. 

Normally, the Earth-occultation data with an elevation an- 
gle below —5° are utilized to inspect the reproducibility of 



Table 3. Signal and background count rates from Suzaku PIN. 



Name 



Count rate (s 'O* 



Sul 



14-20 keV 



30-55 keV 



Observed signal 

NXB 

CXB 
Source net signal 



0.190±0.002 
0.173 
0.011 
0.0055±0.0029 



0.082±0.001 
0.082 
0.002 



Su2 



14-20 keV 



30-55 keV 



Observed signal 

NXB 

CXB 
Source net signal 



0.187±0.002 
0.180 (0.170)t 
0.011 
(0.0068±0.0028)t 



0.081±0.001 
0.084 (0.079)''" 
0.002 



The en'ors are the Icr values, 
t The parenthetic number indicates the corrected value for the 
rescaling factor 



the NXB model. However, Suzaku did not undergo the Earth- 
occultation during our observations. We, instead, compared the 
observed PIN signal in the 30-55 keV range with the modeled 
NXB in the same range, because such PIN data are considered 
to be dominated by the NXB events. Below, we subtracted 
the CXB contribution from the PIN data to derive the correct 
rescaling factor. For Sul, we found no significant difference 
between the data and the NXB model in this energy range 
(the NXB/PIN ratio of 1.5±2.5% in 90% confidence level). 
This consistency allows us to adopt for the PIN data in Sul, 
the NXB model without any correction. However, the mod- 
eled NXB count rate in Su2 is significantly higher than the ob- 
served PIN one by 5.5±2.5%. Such a difference is not reason- 
able from the current uncertainty of the NXB model (typically 
1.3%; Fukazawa et al. 2009). We therefore decreased the NXB 
model count rate by 5.5% for Su2. 

Figure 4 shows the ultimately derived PIN spectra. The 
source spectrum after subtracting NXB and CXB is rep- 
resented in blue. For the spectral analysis, we fo- 
cus on the 14-20 keV band. The net count rates (af- 
ter the NXB and CXB were subtracted) of Sul and 
Su2 are calculated as (5.5±2.9)xl0~'^ counts s~^ and 
(6.8±2.8)xl0"^ counts s"\ respectively (table 3). These 
net count rates are de-convolved to the X-ray flux in 
the same band as (1.2±0.6)x 10^^^ erg cm^^ s^^ and 
(1.4±0.6)xl0^^^ erg cm^^ s^^ for Sul and Su2, respectively. 
Here, the errors of the PIN count rate and the de-convolved PIN 
flux are the la values, in which the statistical and systematic 
uncertainties of the NXB were taken into account. 

Next, we confirm that the majority of the detected (net) 
PIN fluxes are from the two ULXs by inspecting two X-ray 
catalogs; First, we used the Swift BAT 54 months catalog 
(Cusumano et al. 2010), which lists no source in our region, 
and calculated the probability that sources (other than X-1 
and X-2) possibly contributing to the PIN flux are located 
within the PIN view. The detailed procedure is performed 
in the following two steps: (i) The flux limit of the catalog 
is 1x10"^^ erg cm^^ s^^ in 15-150 keV. In order to com- 
pare with the PIN flux limit (^^5x10^^^ erg cm~^ s~^ in 
14-20 keV), we calculated the 14-20 keV flux limit of the 
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Fig. 4. Time-averaged signal spectrum (black) and source (blue) spectrum with Suzaku PIN in (a) Su 1 and (b) Su2. The provided NXB and the simulated 
CXB are represented in red and green, respectively. The NXB for Su2 is corrected for the rescaling factor (§ 4.2.2). The gray shaded range indicates the 
notable energy range of 14-20 keV used for the spectral analysis. 



catalog by assuming a PL shape with F = 2.0. Its value is 
L5xlO~^^ erg cm^^ s^^, which is greater than that of PIN. 
The fact indicates that sources (including the two ULXs) with 
a 14-20 keV flux below the catalog flux limit contribute the de- 
tected PIN flux, (ii) By extrapolating the logN-logS relation of 
the catalog (without Galactic sources at |Z| < 10°), we derived 
that the expected number of sources with a 14-20 keV flux 
between the flux limit of the catalog and that of PIN is about 
2,000 in all sky. From this number and the area ratio of the PIN 
view and all sky, we conclude the probability that such three or 
more sources are simultaneously located within the PIN view 
to be at most 0.1%. 

Second, we used the Chandra catalog (Liu 2011) and ex- 
amined the total 14-20 keV flux of 56 X-ray point-like 
sources (other than X-1 and X-2) within the IC342 re- 
gion. Each source shows the 0.3-8.0 keV flux of 10^^^ to 
10~i5 erg cm~^ s^^. When a PL shape with F = 2.0 is as- 
sumed, the 14-20 keV flux of each source is calculated as 
10~^^ to 10^^^ erg cm~^ s~^. Consequently, its total flux is 
summed up '-^2.9x 10"^'^ erg cm^^ s^^, which corresponds to 
at most 20% of the detected PIN flux. 

4.2.3. XIS+ PIN spectra 

In the XIS analysis (§ 4.2.1), all the XIS spectra were suc- 
cessfully fitted by a phenomenological model (PL, MCD, or 
cutoff-PL). Now, we investigate spectral shapes above 10 keV 
of the ULXs by extrapolating the successful XIS models (with 
the best-fit central parameters in table 2) to the PIN band. We 
note that X- 1 and X-2 are treated simultaneously, since the PIN 
spectrum contains the signals from both ULXs. The PIN effec- 
tive areas to each ULX are much the same"*. We took into ac- 
count the known calibration uncertainty between XIS and PIN 
normalizations"*, namely the PIN normalization is larger than 
the XIS one by a factor of 1 . 164. 

Table 4 shows the spectral model components for X-1 and 

' See "The Suzaku Technical Description" available at 
http://www.astro.isas.jaxa.jp/suzaku/doc/suzaku_td_ao7.pdf for details. 

* See http://www.astro.isas.ac.jp/suzaku/doc/suzakumemo/suzakumemo- 
2008-06.pdf for details. 



X-2. We check the validity of the extrapolation by com- 
paring the modeled flux with the de-convolved PIN flux in 
14-20 keV. For Sul, we examined the PL(X-l)+MCD(X-2) 
combination and the cutoff-PL+MCD one (ii'cut =11 keV 
for X-1 is assumed). Clearly, the PIN spectrum was bet- 
ter reproduced by the former combination (figure 5 a) than 
the latter. This is because the modeled flux of the former 
(0.91 X 10^^^ erg cm^^ s^^) is comparable to the de-convolved 
PIN flux of (1.2±0.6)x 10^12 erg cm^^ s-\ while the flux 
of the latter (0.52x10^^^ erg cm~^ s~^) is smafler than the 
lower limit (Icr) of the de-convolved flux. Consequently, 
the X-1 spectrum in Sul with no turnover (below 20 keV) 
was confirmed as mentioned in § 4.2.1. For Su2, we ex- 
amined the PL and MCD models to X-2 since its XIS 
spectrum was reproduced by both models (§ 4.2.1), while 
only the PL model was considered to X-1. Figure 6 (b) 
clearly indicates that the PL+PL combination with a flux 
of 2.9x10^^^ erg cm^'^ s^^ exceeds the PIN spectrum 
with (1.4±0.7)x 10~^^ erg cm~^ s~^. By replacing the PL 
component with the MCD one for X-2, the model flux of 
1.3x10"^^ erg cm~^ s~^ became consistent with the observed 
one, as displayed in figure 6 (c). 

Based on the present joint analysis, we consider the spectral 
properties as follows: (i) In both observations, the X-1 spectra 
have a PL tail extending up to at least 20 keV. (ii) Both of the 
X-2 spectra cannot be explained by a single PL model. They 
instead have a convex shape. 

4.3. XMM-Newton, Chandra, and Swift Results 

We analyzed the spectra obtained with the three satellites us- 
ing the method same as for the Suzaku XIS (§ 4.2.1). For the 
XMM-Newton EPIC data, we fitted the MOS and pn spectra 
simultaneously. The MOS spectra were multiplied by a con- 
stant factor of 1.0-1.1, because Mateos et al. (2009) reported 
that the MOS flux is higher than pn by 7-9% (<4.5 keV) and 
10-13% (^^4.5 keV). The best-fit parameters are summarized 
in table 5. 

Same as the Suzaku results, most X-1 spectra (XMl, XM3, 
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Fig. 5. XIS plus PIN joint unfolded spectra (pluses). Black and red are respectively XIS spectra of X-1 and X-2, while green is PIN one. (a) The spectral 
models (histograms) are PL for X-I and MCD for X-2 in Sul, (b) PL for X-I and PL for X-2 in Su2, and (c) PL for X-1 and MCD for X-2 in Su2. 



Table 4. Modeled and de-convolved fluxes in 14-20 keV. 



Data 
label 


Model 




Flux in 


14-20 keV 


Notes 


X-1 


X-2 


Model 


NetT 


Sul 


PL 


MCD 


0.91 


1.2±0.6 


figure 5 (a) 


Sul 


cutoff-PL 


MCD 


0.52 


1.2±0.6 




Su2 


PL 


PL 


2.9 


1.4±0.6 


figure 5 (b) 


Su2 

* TU^ riT 


PL 


MCD 


1.3 




1.4±0.6 

,-2 -.-i 


figure 5 (c) 



t The de-convolved PIN flux from the net count rate obtained by subtracting 
the NXB and CXB from the observed PIN count rate. The en'ors are the Icr 
values. 



Chi, Ch2, Swl, and Sw2-|-3) were better reproduced by the PL 
model. They were also reproduced by the cutoff-PL model, but 
we did not derive a significant improvement. Their X-ray lumi- 
nosities (except for Swl) are comparable to that in the Suzaku 
observations ('^5x10"^^ erg s~^). In the Comptonization 
model fitting for these spectra, kT^ and t were coupled, and 
the upper limit of the seed photon temperature was determined 
at 0.5 keV. 

On the other hand, for the two spectra in XM2 and XM4, the 
PL model yielded negative residuals particularly at ^6 keV, in- 
dicating the existence of a turnover (figure 6). In fact, these 
residuals were significantly improved by applying the cutoff- 
PL model (figure 6 and table 5). The reduced-x^(d.o.f.) is from 
1.18(233) by PL to 0.97(232) by cutoff-PL in XM2, and from 
1.12(147) by PL to 0.93(146) by cutoff-PL in XM4, indicating 
the F-test significance of >99.9%. The derived turnover ener- 
gies are £:cut = S.lll J keV in XM2 and S.l^^;^ keV in XM4. 
Their X-ray luminosities ('^1x10*° erg s~^) are higher than 
the other (most) observations with no turnover blow 10 keV 
(or 20 keV with Suzaku) by a factor of 2 to 3. These results 
suggest that the spectral turnover (below 10 keV) is prominent 
only in the higher luminosity phase. In addition, for such spec- 
tra (XM2 and XM4), the Compton parameters were determined 
as kTc ^1.8 keV and r 8.5, although the upper limit of the 
seed photon temperature was 0.5 keV like the other observa- 
tions. 
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Fig. 6. Residuals from the best-fit by the PL model and the cutoff-PL 
model in (a) XM2 and (b) XM4. Black and red represent MOS and pn, 
respectively. 



5. Discussion 



In the present study, we examine the X-ray emission prop- 
erties of the ULX IC 342 X-1, aiming at understanding the na- 
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Table 5. Best-fit parameters for the XMM-Newton, Chandra, and Swift spectra. 



Name 



Data 
label 



Model iVn T or r E cut, kT^, oi kTin 

cm"^) (keV) 



(10' 



,22 , 



0.24^ 

< 0.26 

< 0.02 

-0.02 



Luminosity* 
(10^^ ergs^^^ 



4.0T^ 

3.9 ±0.3 
3.9 ±0.2 



Red-x^ 
(d.o.f.) 



X-1 



XMl 



PL 
cutoff-PL 
compTT 
MCD 



L63±0.09 
1 6+°-^ 



> 7.8 
2.3^ 
1.9±0.1 



1.00(90) 
1.01(89) 
0.97(88) 
1.29(90) 



XM2 



PL 
cutoff-PL 
compTT 
MCD 



0.53 



0.39 ±0.04 

97+0.20 
'-'•^ ' -0.08 

0.15 ±0.02 



2.02 ±0.03 
1 4+01 



8.5 



+0.6 

-1.4 



"•1-1.2 

1.8±0.1 
L54±0.03 



1L1±0.2 
9.4 ±0.4 

Q 9+1.8 
^•^-0.4 

9.2±0.1 



1.18(233) 
0.97(232) 
0.95(231) 
1.72(233) 



XM3 



PL 
cutoff-PL 
compTT 
MCD 



0.30 ±0.03 
0.30 ±0.03 
0.29 ±0.03 
< 0.01 



L82±0.05 

1 0-1+0.05 
J^-»J^-0.08 

5.0"^ 



> 0.01 
4.3^ 
1.78 ±0.05 



4.8±0.1 
4.8±0.1 
5.2±0.1 
4.2±0.1 



1.06(123) 
1.07(122) 
1.08(121) 
2.69(123) 



XM4 



PL 
cutoff-PL 
compTT 
MCD 



0.55 ±0.04 
0.37 ±0.07 
0.46 ±0.05 
0.17 ±0.03 



1.91 ±0.05 
1.2±0.2 
8.4 ±0.9 



1+2.3 
-1.2 
j+0.3 
-0.2 

1.71 ±0.07 



5.V 
1.8^ 



14.2 ±0.4 
12.1±0.2 



1.12(147) 
0.93(146) 
0.93(145) 
1.22(147) 



Chi 



PL 
MCD 
cutoff-PL 
compTT 



0.3 ±0.1 
< 0.05 
0.3±0.1 
0.3±0.1 



1.8±0.2 



1.8 



+0.2 



0.3 

4.5''' 



1.7±0.2 
> 0.01 
5.4^ 



-r 7+0.2 
l-'-O.l 
-I 7+0.2 
1- '-0.7 
4.2+ 



4.7±0.3 
3.7 ±0.2 

^•"-0.7 

4.6±0.3 



0.98(25) 
1.54(25) 
1.02(24) 
1.07(23) 



Ch2 



PL 
cutoff-PL 
compTT 
MCD 



0.3 ±0.1 
3+°-^ 

0.3±0.1 
< 0.08 

U.O_Q 2 

<0.7 
<0.7 
<0.3 



> 0.01 

6.4+ 
1.8 ±0.2 



5.0 ±0.3 
5 0+°-3 
5.0 ±0.3 
4.0 ±0.2 



T2- 

+4^ 
-2 
3 
2 



0.57(27) 
0.59(26) 
0.62(25) 
0.98(27) 



Swl 



PL 
cutoff-PL 
compTT 
MCD 



1.8±0.3 
< 1.9 
11.4+ 



0.001-0.002 
1.4+ 
1.8^ 



"•O-0.3 

<0.9 
<0.8 
<0.3 



}+0.4 

-0.3 



13 



10! 
10 ±1 



-0.7 

0-2 
4+5 

q 7+0.5 
' -0.4 



0.69(29) 
0.64(28) 
0.66(27) 
0.61(29) 



Sw2+3 



PL 
cutoff-PL 
compTT 
MCD 



1.8±0.4 
< 2.2 
4.9+ 



' Same as table 2. 
t Same as table 2. 



> 2.1 
5.3+ 
1 8+°-^ 

-'-■°-0.4 



4.9; 



0.72(12) 
0.78(11) 
0.83(10) 
0.91(12) 



ture of accretion flow in ULXs during the PL state. In con- 
trast with the convex-shaped (thermal) spectra, from which we 
can easily obtain concrete information, such as the mass accre- 
tion rate and the BH spin (as long as the blackbody approxi- 
mation is adopted), it is not so easy to do so for the PL spec- 
tra. Nevertheless, we retrieved some information from the PL 
slope, the existence (or absence) of a turnover, a turnover en- 
ergy (if a turnover exists), its timing properties, and so on. We 
will pay special attention to the turnover energy, since it pro- 
vides a key to understanding the nature of the PL state. 

5.1. Summary of the Observations ofIC342 X-1 

5.1.1. Luminosities and spectral shapes 

We plot the X-ray luminosity (in the energy range of 0.5- 
10 keV) against the observed date and the photon index in 
figure 7. We can see two groups of data occupying sep- 
arate regions, as was reported previously (Feng & Kaaret 
2009; Mak et al. 2011), in these plots: the low luminosity 
group with (4-6) X 1039 erg s"! (Sul, Su2, XMl, XM3, Chi, 
Ch2, and Sw2±3) and the high luminosity group with (1.1- 
1.4)xlO^° erg s"! (XM2, XM4, and Swl). Since both groups 
have been called the "PL" state previously, we expediently 



hereafter call them the low luminosity PL state and the high lu- 
minosity PL state respectively, although the spectra in the high 
luminosity PL state have a turnover (see below). No signifi- 
cant difference was seen in the photon indices between the two 
PL states. The unpublished data (Sul, Su2, Swl, and Sw2±3) 
also follow these trends. Although the strict duty cycle cannot 
be estimated with this sparse data sampling, it seems that the 
low luminosity PL state appears more frequently than the high 
luminosity PL state. Its appearing fraction is low:high = 7:3 
based on the number of the present samples in each state. 

Besides the luminosity differences, X-1 also showed some- 
what different spectral features in hard X-rays between these 
two states. Figure 8 displays the unfolded vFy spectra in the 
unit of keV^ cm~^ s~^ keV~^ in five observations (Sul-2 and 
XM2-4 with good statistics). Evidently, the spectral shapes 
differ at keV. During the high luminosity PL state (XM2 
and XM4), the spectrum has a turnover at around 6 keV. We 
can confirm this property by comparing the spectral fittings to 
the data with the PL model and with the cutoff-PL model. The 
spectra in the high luminosity PL state are significantly im- 
proved, when we apply the cutoff-PL model with the i^-test 
significance at >99.9% (§ 4.3). Whereas, during the low lu- 
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Table 6. Basic properties of the two PL states of IC 342 X- 1 . 



Properties 


Luminosity State 






Low PL 


High PL 


Data set 


Sul, Su2, XMl, XM3, Chi, Ch2, Sw2+3 


XM2, XM4, Swl 


Appearing fraction (present samples) 


70% 


30% 


Range of ix (erg s~^) 


(4-6) X 10^'' 


(l.l-1.4)xl0'"' 


Turnover energy 


Unconstraint (>20 keV) 


~6keV 
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Fig. 7. (a) Long-term X-ray variability, and (b) Lx (X-ray luminosity) versus F (photon index) diagram of IC 342 X-1. The parameters are derived by 
spectral fitting with the PL model. The data taken by the different satellites are indicated by the different colors. 
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Fig. 8. Unfolded spectra of IC 342 X-I in Sul, Su2, XM2, XM3, and 
XM4 shown with different colors. 

minosity PL state (Sul, Su2, and XM3), the spectrum has no 
turnover below 10 keV. In fact, the spectra are not improved by 
the cutoff-PL model (§ 4.2.1 and § 4.3). We estimate that a PL 
tail extends up to at least 20 keV, which we firstly confirmed 
by our Suzaku PIN detector observation (§ 4.2.3). We expect 
similar hard X-ray emissions from this source in the other ob- 
servations during the low luminosity PL state, since the derived 
values (Lx and F in figure 7b) are similar, although they were 
too faint to detect with satellites other than Suzaku. We thus 
conclude that there exists a certain trend in the spectral changes 
with an increase in Lx- Table 6 summarizes the basic proper- 
ties of the two PL states. 



5.7.2. Spectral fitting with the Comptonization model 

The PL-like spectra of Galactic BHBs during the VHS are 
usually explained in terms of the Comptonization model (Done 
et al. 2007). From an overall spectral similarity, we may as- 
sume that the PL-like spectra in both PL states are generated by 
the inverse Compton scattering of soft photons. However, cau- 
tions should be taken here, since there exist some differences 
in the observed properties between the PL states of IC 342 X- 
1 and the VHS of Galactic BHBs; that is, (i) the photon in- 
dices of the PL states {F = 1.6-2.0) are smaller than that of 
the VHS (>2.4; McClintock & Remillard 2006). (ii) The spec- 
tral turnover at a lower energy, around 6 keV, found in the high 
luminosity PL state is not observed during the VHS. (iii) The 
significant variability observed in the VHS is not confirmed in 
the PL states. 

In order to understand the physical condition of the two PL 
states, we show, in figure 9 (a), the best-fit confidence contours 
between r and kF^ derived from the Comptonization model, 
for the five observations (Sul-2 and XM2^). Some specific 
values of Compton y-parameters (y = 0.5, 1, and 2) are also 
indicated in this figure with the dashed lines. 

There are two important facts seen in figure 9 (a): First, 
the two states occupy different regions in this T-kF^ plot: 
fcTc ^1.8 keV and r ~ 8.5 in the high luminosity PL state, 
while kFc (>2 keV) and t (<8) are unconstrained in the low 
luminosity PL state. Second, despite different values of kF^ 
and r, both states show the region with similar y-values; i.e., 
between y = 0.5 and 2. This is a noteworthy feature and its 
reason should be considered carefully, since the y-parameter is 
closely connected to the physical condition of a Comptonizing 
corona (see the next subsection). 
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Fig. 9. (a) Significance contours of tlie best-fit pliysical parameters (t versus kTc) by the Comptonization model for IC 342 X-1. The results in different 
observations are shown with different colors: Sul (black), Su2 (magenta), XM2 (red), XM3 (blue), and XM4 (green). The confidence level is shown for 
the 3<T range. The dashed lines represent the constant Compton {/-parameter of 0.5, 1, and 2. (b) t against kTc of Holmberg XI X-1 by Vierdayanti et 
al. (2010a). The data in black (solid hne) show a generally trend to increase in the optical depth as the X-ray luminosity increases, while the data in gray 
(dotted line) do not. The bolometiic luminosities in black from top to bottom are 11.1, 9.41, 6.92, and 6.63x10^3 erg s"\ and those in gray from top to 
bottom are 5.92, 17.3, 9.27, and ll.lxlO^^ ergs^i. 



Note that PL-like spectra are produced by unsatu- 
rated inverse Compton scattering in the energy range of 
fcTseed hv <^ kTc OR the condition that the seed photon 
temperature (fcTsccd) is much less than the electron temper- 
ature of a Comptonizing corona. This is the case that we 
encounter here, since from the spectral fitting we estimate 
kTscod ^ 0.5 keV, while kT^ > 1 keV in the high luminos- 
ity PL state. A similar argument applies to the low luminosity 
PL state, as well. 

In the following subsection, we interpret the PL states based 
on these spectral properties. 

5.2. Physical Interpretations of the Two PL States 

5.2.1. The high luminosity PL state — possible scenario in- 
voking a Comptonizing corona 

First, we consider the high luminosity PL state, which is 
characterized by a turnover below lOkeV. Kubotaet al. (2002) 
firstly reported the existence of a turnover below 10 keV in 
the PL-shaped spectrum of a ULX (IC342 X-1) with ASCA. 
Miyawaki et al. (2009) also reported a turnover in the PL- 
shaped spectrum of another ULX (M82 X-1) with Suzaku. 
They interpreted that such ULXs are in a state corresponding 
to the "normal" VHS of Galactic BHBs. We, however, already 
pointed out some distinctions between the VHS of Galactic 
BHBs and the high luminosity PL state (see § 5.1.2). 

The PL-shaped spectra with a turnover below 10 keV were 
reported in the high quality spectra of several ULXs analyzed 
by Gladstone et al. (2009). They found that such spectra can 
be explained by the emission from a low-temperature optically- 
thick corona with kTc < 2.5 keV and r > 6. They proposed 
to call it the ultraluminous state (see Soria 2007), claiming that 
the low-temperature optically-thick corona is formed by super- 
critical accretion flows. Their idea has been confirmed with 
global radiation-hydrodynamic simulations by Kawashima et 
al. (2012). Vierdayanti et al. (2010a) also reported the re- 
sults of the PL-shaped spectra with a turnover, which they 



have seen from a long-term monitoring observation of the ULX 
Holmberg IX X-1 (one of the ULXs analyzed by Gladstone et 
al. 2009) taken with XMM-Newton and Swift. They discov- 
ered a general trend of decreasing electron temperatures and 
increasing optical depths as the X-ray luminosity increases, es- 
pecially at T > 6 (the black solid lines in figure 9b). Notably, 
its j/-parameter is ^^1. This r-fcTg relation (including y- 
parameter) is similar to that of IC 342 X- 1 . 

From these similarities, we suggest that 1C342 X-1 is in a 
similar physical condition to that of some ULXs reported by 
Gladstone et al. (2009) and Vierdayanti et al. (2010a). In other 
words, IC 342 X-1 (at least during the high luminosity PL state) 
is likely to undergo supercritical accretion flows, if the interpre- 
tation by Gladstone et al. (2009) and Kawashima et al. (2012) 
is correct. 

However, Vierdayanti et al. (2010a) also reported somewhat 
different behavior First, Holmberg IX X- 1 showed rather con- 
tinuous spectral variations, while IC342 X-1 seems to stay 
in two separate states, although the data sampling is rather 
sparse. Second, spectral changes of Holmberg IX X-1 are 
not always uniquely determined by its X-ray luminosity (see 
the gray dotted lines in figure 9b). It is interesting to point 
that Galactic BHBs exhibit hysteretic spectral transitions; their 
spectral states differ in the rising and decaying phases (see, e.g., 
Miyamoto et al. 1993; Fender et al. 2004). We may expect sim- 
ilar hysteretic behavior, and if so, distinct spectral states appear 
even at similar luminosities. Certainly, we need more continu- 
ous observations to see if these differences really exist. 

We briefly comment on the case of the Galactic BHB 
GRS 1915+105, in which the mass of the BH is known to be 
14±4 Mq (Greiner et al. 2001). At near the Eddington lu- 
minosity (apparent L ^ 0.7iEdd) under the supercritical ac- 
cretion situation, the source showed a unique state with a low 
electron temperature {kTc ^ 3 keV) and optically-thick (r ^ 9) 
Comptonizing corona (Vierdayanti et al. 2010b). Such features 
are quite reminiscent of those of the PL state with a turnover of 
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ULXs, and this similarity may suggest similar physical condi- 
tions; i.e., the supercritical accretion may also be occurring in 
the ULXs. We need to point an annoying fact, however; that is 
the seed photon temperatures are distinct: GRS 1915+105 ex- 
hibits a disk-dominated spectrum with kTgccd ^ 2 keV, while 
our data of IC 342 X- 1 show much lower seed photon temper- 
atures, below 0.5 keV. The origin of this discrepancy is left as 
an open question. 

5.2.2. A possible theoretical interpretation of the high lumi- 
nosity PL state 

Here, we try to understand a physical cause of yielding 
y ~ 1 observed in the high luminosity PL state based on 
the framework of the disk-corona model. Let us consider 
the energy balance of the corona and the accretion disk. 
The former can be described by (3+,o„a = Q^rona' where 
Qcorona Qcorona represent the heating and the cooling 
rates of the corona, respectively. The latter is approximated 
by Qmrona ^ {y/2)cUsccd, whcrc C/sccd = aT^ccd is the seed 
photon energy density generated on the surface of the disk, 
and a and c are the radiation constant and the speed of light, 
respectively (see, e.g., Meyer et al. 2000). From the en- 
ergy balance within the disk, on the other hand, we have 
Qisk = Qdisk = '^T^ccd = c?/sccd/4, where Q+^^, Q-.^^, and a 
are the heating rate, the cooling rate of the disk, and the Stefan- 
Boltzmann constant. Suppose that most accretion energy liber- 
ated within the disk is transported to the corona by the magnetic 
field (e.g., Haardt & Maraschi 1991). Since about a half of the 
radiation scattered in the corona heats up the disk material, the 
energy balance in the disk leads to QtJorona/^ ^ Qdisk- Finally 
we obtain y ^ 1. Note that the above calculations can apply 
only when a Comptonizing corona is optically thick (r ^ 1). 

5.2.3. The high luminosity PL state — possible scenario in- 
voking a disk 

In the above scenario, the existence of a low-temperature 
optically-thick corona was assumed. Inversely, we here con- 
sider another possibility; that is, the spectral turnover in the 
high luminosity PL state may be given by direct disk emissions. 
The origin of the turnover cannot be the standard disk, since 
both spectra in XM2 and XM4 (namely the high luminosity PL 
state) were unaccountable by the MCD model (table 5). Hence, 
we examined the slim disk emission by applying an extended 
MCD model (so-called the p-free model) with the additional 
parameter p, which assumes a flatness of the temperature (T) 
profile of the disk; i.e., T{r) oc r^P, where r is the radius of the 
disk. The extended MCD model with p = 0.75 is equal to the 
"normal" MCD model. If the spectra are from the slim disk, 
it is expected that p approaches from 0.75 to 0.5 as the X-ray 
luminosity increases (e.g., Watarai et al. 2001). The same con- 
clusion was derived, even if mass loss by a radiation-pressure 
driven outflow is taken into account (Takeuchi et al. 2009). In 
fact, such a p-Lx variation was reported in a BHB (Kubota & 
Makishima 2004) and a ULX (Isobe et al. 2009) in the slim disk 
state. For IC 342 X-1, the individual spectra can be reproduced 
by the extended MCD model (reduced-x^ = 0.95 in XM2 and 
0.92 in XM4), but the derived p-Lx variation (p = 0.54±0.01 
and Lx = 1.0x10"^" erg s^^ in XM2 to p = 0.57±0.02 and 
Lx = 1.2x10'"' erg s^^ in XM4) does not agree with the ex- 
pected one. We thus conclude that, at least, the 1C342 X-1 
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Fig. 10. Spectral energy distributions for tlie Comptonization model 
with parameters of each PL state. The black solid line shows the high 
luminosity PL state with kTc = 1.8 keV and r = 8.5. The gray lines 
show the low luminosity PL state with kTc = 1 keV and r = 3.6 (solid) 
and kTc = 20 keV and r = 1.8 (dashed). The seed photon temperature 
of both states is assumed to be 0. 1 keV. The range used in the present 
paper is shown with gray shading (0.5-10 keV and 14-20 keV). The 
energy range of the Astro-H and the NuSTAR satellites is shown at 
top-right. 

spectra in the high luminosity PL state are from neither the 
standard disk nor the slim disk. 

5.2.4. The low luminosity PL state 

Next, we discuss the low luminosity PL state, which is char- 
acterized by a PL tail extending up to at least 20 keV. Since 
the low luminosity PL state also shows y ^ \, like the high 
luminosity PL state, these two PL states can be described by 
a similar corona model with similar energy balance. In other 
words, the distinction between the two PL states can be under- 
stood only by different coronal temperatures or coronal optical 
depth. 

Two unresolved issues remain: First, it is not clear which of 
kTc or T is the principle (controlling) parameter solely from the 
spectral fitting results. If we stand on a scenario with the su- 
percritical accretion rate, increase in t is a cause of the electron 
temperature, since the amount of outflow gas (i.e., t) can vary, 
depending on the disk luminosity (Kawashima et al. 2009). 
The observed behavior of the two PL states is consistent with 
their calculation. 

Second, the electron temperatures of the Comptonizing 
corona (or outflow, if it every exists) in the low luminosity PL 
state has not yet been precisely measured so far. If the low lu- 
minosity PL state is somewhat similar to that of the VHS of 
Galactic BHBs, we expect higher electron temperatures, say, 
kTc ^ 20 keV. Or alternatively, there may be no states like 
the VHS in ULXs, which is supported by the absence of large 
luminosity variations in ULXs. In order to settle down this is- 
sue, more sensitive, higher energy (>20 keV) observations are 
required. We illustrated the expected spectra of two represen- 
tative cases in figure 10. In one case, there is a Comptonizing 
corona (or outflow) with a similar coronal electron temperature 
to that of the VHS (fcTc = 20 keV), and hence the turnover en- 
ergy is ~60 keV. While in the other case, the electron tempera- 
ture is much less, kT^. = 7 keV, leading to a lower turnover en- 
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ergy is '^20 keV. Here, we fixed the optical depths to be r = 1.8 
and 3.6, respectively. Such observations will be made possible 
by the launch of Astro-H (Kokubun et al. 2008; Takahashi et 
al. 2010) and NuSTAR (Harrison et al. 2005). 

5.2.5. Relation with the convex state 

Finally, we briefly mention the relation between the PL 
states and the convex state. In 1C342 X-1, the convex state 
has been observed on 1993 September with ASCA (Okada et 
al. 1998; Kubota et al. 2001). Its spectral shape is apparently 
different from those in the PL states. In fact, the spectrum in 
the convex state can be reproduced by the MCD model, which 
does not reproduce the spectra in the PL states. However, the 
0.5-10 keV luminosity in the convex state (^1.3 x IC*" erg s^^ 
assuming 3.3 Mpc) is comparable to that in the high luminosity 
PL state. That is, two distinct states appear at a similar lumi- 
nosity. These facts imply that the spectral state of ULXs (at 
least IC 342 X- 1 ) is not uniquely determined by only the lumi- 
nosity (or the accretion rate) but is determined in a more com- 
plicated fashion. One possibility is a hysteretic relation known 
in Galactic BHBs, as we mentioned earlier (§ 5.2.1). More fre- 
quent observations are required to elucidate complex spectral 
variations. 

6. Summary 

We have analyzed X-ray spectra of the unique ULX IC 342 
X-1 by using our two Suzaku data and archival data from mul- 
tiple XMM-Newton, Chandra, and Swift observations. X-1 
clearly showed two distinctive sub-states: (i) The low lumi- 
nosity PL state appears in 70% of all the present observations. 
The X-ray luminosities are (4-6) xlQ-^^ erg s^^ (0.5-10 keV), 
and the spectra have a PL tail extending up to at least 20 keV. 
The coronal parameters derived from a Comptonization model 
are strongly coupled {kT^ > 2 keV and r < 8). (ii) The high 
luminosity PL state appears in 30%. The X-ray luminosities 
are (l.l-1.4)xlO''° erg s~^ (0.5-10 keV), and the spectra have 
a turnover at about 6 keV. The coronal parameters are deter- 
mined at fcTo 1.8 keV and r ~ 8.5. 

Since the two PL states showed the similar Compton y- 
parameter of y ~ 1, we suggest that the two PL states are de- 
scribed by a similar corona model with similar energy balance. 
The electron temperature and the optical depth respectively de- 
creases and increases as an increasing X-ray luminosity. This 
behavior of IC 342 X-1 is similar to that of some ULXs (espe- 
cially HolmberglX X-1) probably under the supercritical ac- 
cretion situation, or that of the Galactic BHB GRS 1915+105 
in the unique state with the supercritical accretion. These facts 
imply a possibility that the two PL states of IC342 X-1 are 
formed by the supercritical accretion flows. 

On the other hand, an unresolved issue for the relation be- 
tween the high and low luminosity PL states remains. More 
sensitive and higher energy observations achieved by the next 
generation X-ray satellites Astro-H and NuSTAR would unveil 
the detailed physical condition of the low luminosity PL state. 
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